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We consider the inverse problem of in-line holography, applied to minimally-destructive imaging 
of cold atom clouds. Absorption imaging near resonance provides a simple, but destructive measure- 
ment of atom column density. Imaging off resonance greatly reduces heating, and sequential images 
may be taken. Under the conditions required for off-resonant imaging, the generally-intractable in- 
verse problem may be linearized. A minimally-destructive, quantitative and high-resolution image 
of the atom cloud column density is then retrieved from a single diffraction pattern. 

PACS numbers: 42.30.Rx, 42.40.Ht, 32.80.Pj 
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The simplest optical measurement of object structure 
is made by illuminating the object with radiation and 
recording the diffraction pattern produced (Fig. The 
inverse problem of retrieving the structure of a non- 
crystalline object from its Fresnel diffraction pattern 
has been studied since Gabor's first incomplete solution, 
known as in-line holography P]. This Communication 
solves the particular inverse problem of retrieving the 
structure of a cold atom cloud from a single diffraction 
pattern. 

Minimally-destructive imaging of cold atoms requires 
weak absorption. High-resolution imaging also places 
constraints on phase shifts. Under these assumptions, we 
derive a linear solution which retrieves the column den- 
sity of the atom cloud from a single diffraction pattern. 
We apply this solution to demonstrate off-resonant imag- 
ing of a cold atom cloud without beamsplitters, phase- 
plates, or imaging optics. 

Gabor's in-line holography recovers an approximation 
of the original wavefield by illuminating a photographic 
transparency of the diffraction pattern. The recon- 
structed wavefield is contaminated by the superimposi- 
tion of an out-of- focus twin image Other forms of 
holography use a reference beam to record an interfer- 
ence pattern, rather than a diffraction pattern, and so 
separate the twin image 0. In this Communication we 
demonstrate a non-holographic method of retrieval. Such 
methods have been proposed when it is inconvenient or 
impossible to generate a coherent reference beam. The 
first step is common with Gabor's method: a diffrac- 
tion pattern is recorded without the need for optics such 
as lenses and beamsplitters. In the second step, rather 
than reconstructing the wavefield by optical propagation 
(physical or numerical) , an image is extracted by solving 
an inverse problem with specified constraints. 

There is insufficient information in a single intensity 
image to retrieve both the amplitude and the phase of 
the wavefield. This information deficit may be balanced, 
and the inverse problem solved, if the object is assumed 
to be purely absorbing ^ or purely phase-shifting 0, 
but these assumptions are seldom valid in practice. 

Instead, we present a single-image solution based on 



the assumption of a monomorphous object (one made of 
a single material), so that both the phase-shift (p and 
the absorption fi of the object are proportional to the 
column density of material along the optical path p(x) = 



/i(x) — kj3p('x) and 0(x) ~ kSp{x.). 



(1) 



The variable x represents coordinates in a plane trans- 
verse to the incident wave propagating along the z-axis 
and k — 27r/A is the wavenumber for illuminating radi- 
ation of wavelength A. The absorption and phase coef- 
ficients P and S correspond to a refractive index of the 
form 



n = l + N{r){d + il3), 



(2) 



with N(r) the atom number density. This mono- 
morphous object assumption has been used in compen- 
sating defocus and spherical aberration in electron mi- 
crographs and transport-of-intensity imaging ^ . 

Immediately after an optically thin object (Fig. Q^a)), 
an incident scalar plane-wave of amplitude /o becomes 
/(x) = /o exp ( — /i(x)+i(/)(x)) , and this wavefield may be 
propagated through a distance z by the Fresnel transform 
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FIG. 1: Recording a diffraction pattern. Vertical lines repre- 
sent intensity profiles. In part (b), point-source illumination 
magnifies the diffraction pattern. 
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/(x)exp(^lx-x'p) dx'. 



(3) 

The Fresnel approximation agrees closely with the com- 
plete scalar diffraction theory, except for propagation at 
large angles to the axis or within a few wavelengths of 
the object. Optical detectors measure intensity / = |/p 
and it can be shown that the Fourier transform JFof 
the diffracted intensity measured at z can be expressed 
in terms of the object-plane wavefield /(x) as 

/+00 
r(x+Azu/2)/(x-Azu/2) 
-oc- 

X exp(— 27rjx • u) dx, (4) 

in which u is the spatial frequency conjugate to x. Writ- 
ten in terms of absorption and phase-shifts, this is 

/+00 
exp{— /i(x+ Azu/2) — /i(x — Azu/2) 
-oc- 

+i [^(x-Azu/2) - 0(x+Azu/2)]} 
X exp(— 27rjx • u) dx. (5) 

Assuming both real and imaginary parts of the expo- 
nential are small, we expand and apply the Fourier shift 
theorem to yield: 

T[I{x,z)] = /o((5(u)-2cos(7rAzM^)J^^(x)] 

+2sin(7rAzu2)jC|0(x)]). (6) 

This expression |^ relates absorption and phase-shift 
to the intensity of the diffraction pattern. The linearizing 
assumption used in obtaining Eq. © implies: 

2Ai(x) < 1 (7) 
and |0(x-t-Azu/2) - (/«(x- Azu/2)| < 1. (8) 

The object must not be strongly absorbing, but it need 
not be completely transparent. The phase-shift should 
obey the finite difference condition Eq. (jHJ , which re- 
stricts large variations in the phase-shift to coarse struc- 
tures in the object. Note that weak phase-shift (|(/'(x)| ^ 
1) is sufficient to satisfy Eq. ((HJ but is not necessary 
This phase condition may always be met at small z, but 
phase objects of many radians thickness may require im- 
practically small propagation distances, and phase shifts 
of order 1 radian are preferable. 

For monomorphous objects obeying Eq. there is 
then a linear shift-invariant relation between the normal- 
ized contrast I / Iq — 1 and the column density p: 



I-Iq 
lo 



2k ((5 sin(7rAzu^ 



(9) 




FIG. 2: The contrast transfer function h, normalized, for 
phase-advancing (solid line) and phase-retarding (dashed line) 
monomorphous objects. At low spatial frequencies the con- 
trast approaches the value for an in-focus absorption image. 



The factor h(u: z) = (5sin(7rAzM^) — (3 cos{tt\zu^) is 
termed the contrast transfer function (CTF), and is plot- 
ted in Fig.|21for positive and negative values of the phase 
coefficient S. Equation IHl can be solved formally for p, 
but the zeros in the CTF render the retrieval an ill-posed 
inverse problem. 

The inverse problem may be regularized, for exam- 
ple by the Tikhonov method Rather than dividing 
J-[I/Iq — 1] by h{u;z), the Tikhonov filter retrieves the 
column density by the modified division 



P(x) 



2k 



h{u; z) 



h'^{u; z) + 
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I-Iq 
lo 



(10) 



which closely approximates division by the CTF except 
at spatial frequencies where the CTF is near zero. Larger 
values of the Tikhonov parameter a reduce the amplifica- 
tion of noise in the retrieval process, but at the expense 
of image distortions. Smaller values yield less distorted 
but noisier retrievals; a normalized alpha value of 0.2 was 
used. Algorithmic optimization of a is possible, for exam- 
ple using Fourier- wavelet regularized deconvolution 0| . 

It is clear from Fig. |21 that the solution is more sta- 
ble if the object advances the phase of the incident wave 
{S < 0), and the CTF zero-crossing at low spatial fre- 
quencies is avoided. The column density may be re- 
trieved for phase-retarding objects {S > 0) but, as shown 
in Fig. [3 the focusing action of the phase-shift cancels 
the absorption contrast at low spatial frequencies and 
lower quality retrievals result. 

If the diffraction pattern is re-imaged by a lens, the 
system may be defocused behind the object so that the 
effective propagation distance z is negative. It follows 
from Eq. (O that the sign-condition on S is then reversed. 
For lensless imaging, negative z cannot be achieved and 
the object should be phase-advancing. 

Magnified images can be retrieved even without lenses. 
Rather than using plane-wave illumination, a point- 
source of light a distance Ri before the object produces 
a spherical wave incident on the object (Fig. ^b)). The 
diffraction at detector distance R2 is magnified by the 
geometric factor M — {Ri + R2)/Ri, but is otherwise 



3 



identical to the plane-wave pattern of Fig. HTa) at the 
effective propagation distance ZeS = R2/M 

Conventional optical materials are phase-retarding 
[5 > 0) but for x-ray imaging [3, and for imaging atomic 
gases with light blue-detuned from an atomic resonance, 
the phase is advanced. We now show that our solution to 
the diffraction imaging inverse problem is exactly suited 
to off-resonant imaging of cold atom clouds. 

To date, all measurements of ground-state BECs have 
been made with near-resonant optical probes. On- 
resonance absorption imaging is destructive for most 
BEC configurations. Imaging with an off-resonant probe 
reduces heating due to spontaneous emission, with the 
cloud instead shifting the phase of the probe beam. Dark- 
field Q and Zernike phase-contrast techniques of 
phase microscopy have been used to render these phase- 
shifts visible and hence obtain sequential, minimally- 
destructive images of BEC. Other methods investigated 
for off-resonant imaging include an intcrfcrometric tech- 
nique equivalent to off-axis image holography and 
a propagation-contrast method based on transport-of- 
intensity |T^ . Although minimally-destructive imaging 
has been crucial to observing many dynamic processes in 
BEC, destructive absorption imaging is still widely used. 

As shown above, free space propagation produces 
phase contrast without optics. Gaussian fits to atom 
cloud images (peak column density and diameters only) 
have been extracted from fitting diffraction patterns [13 ■ 
Our solution to the inverse problem retrieves detailed 
column-density images of cold atom clouds, without re- 
quiring Zernike phase-plates or interferometry. 

In the two-level approximation, the refractive index of 
an atomic gas is 



l + iV(r) 



ctqA i - 2A 



(11) 



where N is the number density of atoms, A is the 
detuning in natural linewidths and ctq is the resonant 
cross-section (3A^/27r for closed transitions). Compar- 
ison with Eq. Q confirms that such an atomic cloud 
is a monomorphous object, with absorption and phase- 
shifts through the cloud given by Eq. QJ. Provided 
that the atom cloud satisfies the weak-absorption con- 
dition Eq. O and limited-phase condition Eq. ||SJ), the 
CTF relation Eq. ^ applies. Substituting the /3 and 6 
coefficients from Eq. yields the cold atom CTF 



h{u; z) 



2fc(I + 4A2 



■ [2A sin(7rAzw^) — cos(7rAzu^)] , 



(12) 

which is then used in Eq. (|I0|) to retrieve the column den- 
sity of the atom cloud. If the detuning is blue of reso- 
nance, the atom cloud advances the phase of the incident 
light, and the low-frequency CTF zero is avoided. 

An important feature of this application to cold atom 
imaging is the regularizing effect of residual absorption. 



At zero spatial frequency, the CTF falls to ^ 1/2 A of 
its maximum value, but does not vanish completely as 
it would for a pure phase object. Even small residual 
absorption is sufficient to stabilize the retrieval and then 
the Tikhonov modified form Eq. H1U|I need only be used 
at higher spatial frequencies above Umin = 1/ V^Xz. Such 
partial regularization greatly reduces distortion, and re- 
trievals approach the optimal linear estimate (Wiener 
filter) which can only be calculated with full advance 
knowledge of the object power spectrum. 

This linearization of the inverse problem is only valid 
if the atom cloud meets the absorption and phase con- 
ditions H7I8|) . Minimally-destructive imaging necessarily 
obeys the weak-absorption condition. The phase con- 
dition Eq. JSJ also broadly concurs with physical con- 
straints due to refraction and resolution. Light refracted 
by the object must remain within the numerical aperture 
of the lens (or detector) Q . It follows that structures at 
the diffraction limit of the imaging system should have 
phase variations less than one radian. Objects satisfying 
this 'thin object' condition - that the detailed structure 
of the object must vary by less than a radian ~ are likely 
to also satisfy the slowly- varying phase condition Eq. |(SJ • 

Detuning the probe light by /9maxCo/4 full linewidths 
from resonance reduces the peak phase-shift to order one 
radian. BECs typically have resonant optical-densities 
PmaxCo ~ 300 and so detunings must be of order lOOF 
to meet the refraction condition. At such detunings, 
many images may be taken before the cloud is appre- 
ciably heated. In the shot-noise limit, further increasing 
detuning and intensity does not improve the SNR beyond 
a limiting value, and in the presence of technical noise will 
reduce the SNR. It has been shown that this SNR limit is 
determined only by the number of spontaneous emission 
events and condensate parameters [l9j- 

In proof-of-principle experiments, the point-projection 
configuration of Fig.^b) was used to image a cold atom 
cloud of peak resonant optical density 2.2. A weak 
linearly-polarized probe beam, detuned -1-3. IF from the 
^^Rb 4Si/2(F = 3) 5P3/2(i^ = 4) transition, diverged 
from the cleaved endface of a single-mode optical fiber, 
expanding for i?i = 125 mm before passing though a 
cloud of cold atoms around 300 /im in diameter held in 
a magneto-optical trap (MOT) vacuum chamber. Trap- 
ping beams were turned off for the 300 fj,s duration of the 
probe pulse. The beam propagated a further distance 
i?2 = 155 mm to a CCD camera (Photometries Cool- 
snap HQ, 1392x1040 pixels, 6.45 /zm pitch), producing 
the diffraction pattern shown in Fig. |31 The column- 
density image was retrieved using Eq. 1)10(1 ; retrievals take 
around one second on a Pentium-Ill processor using stan- 
dard discrete Fourier transform algorithms [2ll |. 

It is usually not possible to measure the propagation 
distances i?i and R2 accurately enough to produce an op- 
timal retrieval. Instead, the retrieval is performed with 
the contrast transfer function h(u; z) evaluated at vari- 
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FIG. 3: Above, the recorded diffraction pattern. Below, the 
column density retrieved using Eq. IIOII . shown at twice the 
magnification of the diffraction pattern. Plots averaged over 
the central five rows of pixels are shown below the images. 
The overplotted gray trace shows the column density calcu- 
lated from an on- resonance in-focus absorption image. 



depth-of-field and post-hoc focusing are retained. 

In practice, the propagation distance is constrained by 
resolution limits. For a detector of diameter D, the min- 
imum resolvable line spacing in the retrieved image is of 
order 2XR2/D, as it is for a lens of the same diameter 
in the same position. The Ri distance should then be 
chosen to provide sufhcient magnification that resolution 
is limited by diffraction and not by the pixel size. In 
our experiment, optical access limited the resolution to 
30 ^m. We can predict that a BEC in a glass cell imaged 
with i?i = 12 mm and R2 — 60 mm on a D = 25 mm 
CCD with 9 /im pixels yields a pixel resolution of 3 /im 
and a diffraction-limited resolution of 3.7 /im. Further, 
the CTF depends on z, A, A and u. Once z is set by 
'focusing', the remaining quantities are readily measured 
to better than 1%. Such well-defined parameters and the 
lack of lens aberrations yield highly quantitative column 
density measurements. 

We solve the inverse problem of retrieving a quan- 
titative column density image from a single diffraction 
pattern by exploiting the proportionality between ab- 
sorption and phase shift through a single-material ob- 
ject. The predicates of the solution are uniquely suited 
to imaging cold atom clouds. Lens aberrations are pre- 
cluded by avoiding the need for image formation entirely. 
Beam-splitters, phase-plates and other optical elements 
are also unnecessary. The holographic record allows re- 
focusing after the image has been acquired. The solution 
can be used with existing absorption imaging systems 
simply by defocusing the imaging lens. We calculate 
near- wavelength resolution when using point-projection 
to image Bose-Einstein condensates. 



ous values of z until a sharp image is retrieved. Thus 
focusing is performed in software when retrieving, rather 
than by adjusting lens positions when imaging. As a 
corollary, the retrievals show holographic depth-of-field: 
one diffraction pattern can be used to retrieve images at 
many different z values. The very real problems of focus- 
ing the optics, and of the limited depth-of-ficld inherent 
in high-resolution imaging, are completely obviated. 

While the precise setting and knowledge of the propa- 
gation distance is immaterial, its coarse setting affects the 
shape and contrast of the diffraction pattern and hence 
the signal-to-noise ratio (SNR) of the retrieved image. At 
short distances only residual absorption contrast will be 
rendered. At large propagation distances and for small 
phase-shifts, the root-mean-square SNR approaches 71% 
of that obtained with the Zernike technique JQj . 

A further advantage of the point-projection configura- 
tion is the absence of lenses, and their resolution-limiting 
aberrations. Diffraction contrast may also be produced 
by defocusing an existing absorption imaging system, 
which may be more convenient than placing the camera 
very close to the object. The advantages of holographic 
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